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ABSTRACT 

We estimate minimum energy magnetic fields (B m [ n ) for a sample of galaxies with measured gas surface 
densities. The sample spans more than four orders of magnitude in surface density from normal spirals to 
luminous starbursts. We show that the ratio of the minimum energy magnetic pressure to the total pressure in the 
ISM decreases substantially with increasing surface density. For the ultra-luminous infrared galaxy Arp 220, 
this ratio is ~ 10~ 4 . Therefore, if the minimum energy estimate is applicable, magnetic fields in starbursts are 
dynamically weak compared to gravity, in contrast to our Galaxy and other normal star-forming spiral galaxies. 
We argue, however, that rapid cooling of relativistic electrons in starbursts invalidates the minimum energy 
estimate. We critically assess a number of independent constraints on the magnetic field strength in starburst 
galaxies. In particular, we argue that the existence of the FIR-radio correlation implies that the synchrotron 
cooling timescale for cosmic ray electrons is much shorter than their escape time from the galactic disk; this in 
turn implies that the true magnetic field in starbursts is significantly larger than B m [ n . The strongest argument 
against such large fields is that one might expect starbursts to have steep radio spectra indicative of strong 
synchrotron cooling, which is not observed. We show, however, that ionization and bremsstrahlung losses 
can flatten the nonthermal spectra of starburst galaxies even in the presence of rapid cooling, providing much 
better agreement with observed spectra. We further demonstrate that ionization and bremsstrahlung losses are 
likely to be important in shaping the radio spectra of most starbursts at GHz frequencies, thereby preserving 
the linearity of the FIR-radio correlation. We thus conclude that magnetic fields in starbursts are significantly 
larger than B m ; n . We highlight several observations that can test this conclusion. 

Subject headings: galaxies:general — galaxies:magnetic fields — galaxies:starburst — infrared:galaxies — 
radio continuum:galaxies 



1. INTRODUCTION 

The magnetic energy density of the Galaxy is observed to 
be in rough equipartition with the cosmic ray energy density 
and the turbulent pressure. The sum of these yields a total 
midplane pressure consistent with that required by hydrostatic 
equilibrium, given the mass distribution in the solar neighbor- 
hood (Boulares & Cox 1990). This equipartition magnetic 
field is roughly 6 /xG at the solar circle (e.g., Beck 2001). 
Similar magnetic field strengths are found in most normal 
star-forming spiral galaxies (Fitt & Alexander 1993; Niklas 
1995). 

Starbursts and ultra-luminous infrared galaxies (ULIRGs) 
have much higher surface densities and turbulent velocities 
than local spirals. As an example, Arp 220 has a gas sur- 
face density ~ 10 4 times higher than the Galaxy (Downes 
& Solomon 1998). If Arp 220 has a magnetic energy den- 
sity in rough equipartition with its hydrostatic pressure, the 
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implied magnetic field strength would be ~ 0.03 G on few- 
hundred parsec scales (Fig.Q Table|2j. This is substantially 
larger than the value of ~ 1 mG typically inferred in ULIRGs 
(e.g., Condon et al. 1991) using the observed radio emission 
and the classic "minimum energy" argument (Burbidge 1956; 
Longair 1994). 

Throughout this paper we use the term "minimum energy" 
to refer to the magnetic field strength inferred using the ob- 
served radio flux and assuming comparable cosmic ray and 
magnetic energy densities (B m i n ; eq. Q). We reserve the term 
"equipartition" for magnetic field strengths that are dynami- 
cally important with respect to gravity (B eq ; eq. (3)). 

We find that the large discrepancy between the minimum 
energy and equipartition magnetic field estimates in Arp 220 
is generic to starburst galaxies; this discrepancy motivates 
the analysis of magnetic field strengths in luminous starbursts 
presented in this paper. In 5j2]we present a sample of galax- 
ies with measured gas surface densities and radio fluxes, for 
which both minimum energy and equipartition field strengths 
can be estimated. In |3]we summarize arguments in favor of 
B ~ B min and against significantly larger fields. Rebuttals to 
a subset of these arguments are given in §0] where we also 
provide independent arguments why magnetic fields in star- 
bursts are likely to be 3> B m [ n . These arguments draw heavily 
on the FIR-radio correlation for star-forming galaxies and the 
observed radio spectra of starbursts. Finally, in |5]we sum- 
marize our results and discuss their implications. 
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FlG. 1. — Minimum energy magnetic field (B mm ', eq. QJ) versus measured gas surface density (T, g ). Normal star-forming galaxies (open squares; Table 0, 
starburst galaxies (filled squares; TablelH. and the Galaxy (B m j n = 6 /^G at the solar circle, e.g., Beck 2001; S g 



: 2.5 x 10~ 3 g cm" 2 



Boulares & Cox 1990) are 



shown. To indicate the uncertainties at high T} g , the open triangles show B mm for IC 883 (Arp 193), Mrk 273, and the individual nuclei of Arp 220, as inferred 
using Eg and the radial size from Downes & Solomon (1998) (their "extreme" starbursts; see their Table 12) and with the radio flux from Condon et al. (1991) at 

8.44 GHz (see Table |3*1. The dotted line is the equipartition magnetic field (eq. |3)) and the dashed line is the scaling B m j n oc Y? g (eq. II 11 1 The latter scaling 
is expected if the cosmic ray electron cooling timescale is much shorter than the escape timescale from the galactic disk, in which case the true magnetic fi eld is 
also significantly larger than the minimum energy estimate (see eq. 1111 ). The solid lines show the scalings B oc Y," (with a = 0.9, 0.8, and 0.7) used in 614.31 



2. DATA & RESULTS 

Longair (1994) derives the minimum energy magnetic field 
strength as (Vol. II, eq. [19.30]) 



(1) 



J 2 ^v,23 "6 "GHz 

where L„ i23 = L„/10 23 W Hz" 1 , ^ GHz = v/10 9 Hz, V 6 = 
V/(100) 3 pc 3 is the volume of the emitting region, <52 = <5/10 2 
is the ratio of the energy in cosmic ray ions to the energy in 
cosmic ray electrons, and the derivation assumes that the cos- 
mic ray electron energy spectrum is n(pf) oc 7~ p with p = 2.5. 
We have scaled the results above for parameters appropriate 
to compact starbursts. For parameters typical of the Galaxy, 
~ 3 x 1(T 2 and V 6 - 7 x 10 4 so that B min w 10" 5 G. 

The minimum energy magnetic field strength is often re- 
ferred to as "equipartition" because the corresponding mag- 
netic energy density (UB min ) is approximately equal to the total 
cosmic ray (electron + ion) energy density if B m \ n actually 
obtains. In this paper, however, we use the term "equiparti- 
tion" to denote magnetic energy densities comparable to the 
total hydrostatic pressure of the ISM. For a gas disk of sur- 



face density S g , hydrostatic equilibrium with the self-gravity 
of the disk implies a midplane pressure 

^2 



(2) 

We measure equipartition via the parameter r\ defined by Ub = 
■qP. The equipartition (77 = 1) field strength is thus 



B eq «(8tt 2 G) 1/2 E„ «2E„ mG, 



(3) 



where T, g is in cgs units. Theoretically, there are significant 
uncertainties about the origin of magnetic fields in galactic 
disks, and so it is difficult to estimate the expected value of ij. 
Field strengths as large as equipartition are, however, possi- 
ble if the magnetic energy density equilibrates with the turbu- 
lent energy density of the ISM. The scale-height h of a self- 
gravitating galactic disk is given by h ~ 5v 2 /2nGYi g , where 
£v is the turbulent velocity. As a result, the energy density in 
turbulent motions is pSv 2 ~ 7rGE 2 , where p is the mean den- 
sity of the ISM. Thus, if the field is amplified by the turbulent 
motions so that B 2 /8tt ~ p8v 2 , equipartition follows. 

The gravitational potential of normal star-forming galax- 
ies includes a significant contribution from stars (and per- 
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haps dark matter) on the scale of their effective radii. Ne- 
glecting for simplicity the different scale-heights of stars and 
gas, equation Q should then be replaced by P « 7TGEgE tot w 
TrG^f~ l , where f g is the gas fraction and E tot is the total 
surface density. Typically, f g is in the range ~ 0.1 -0.2 for 
normal galaxies and f g ~ 0.5 for starbursts. Thus, the mag- 
netic field could in principle be larger than equation (|3j by a 

-1 /2 

factor of f g ' and still be sub-dominant with respect to grav- 
ity. In this paper, these distinctions are not important because 
we focus on gas-rich starbursts and use equation l|3} primar- 
ily as a useful point of comparison for field strengths inferred 
using the minimum energy argument (eq. 0). 

Figure shows the minimum energy magnetic field Z? m i n 
inferred using equation Q as a function of the surface density 
E g of the galactic disk for a sample of galaxies ranging from 
local spirals to ULIRGs. The sample is based on galaxies 
with surface densities and sizes used by Kennicutt (1998) to 
study the global properties of star formation across a range of 
galaxies (the "Schmidt Laws"). The properties of the normal 
star-forming and starburst galaxies are listed in Tables \l\ and 
12 respectively. Systems described as starbursts by Kennicutt 
are denoted here by solid squares, whereas his normal star- 
forming galaxies are labeled by open squares. 

In computing the radio emitting volume (eq. Q), we as- 
sume that the scale height of the synchrotron emission is 
hrad = 500 pc for the normal galaxies (Beuermann et al. 1985; 
Dumke & Krause 1998) and h ra d = lOOpc for the starbursts 
(filled squares) (e.g., Klein et al. 1988; Seaquist & Odegard 
1991; Condon et al. 1991). 8 We use radial sizes from Kenni- 
cutt (1998); for normal galaxies, these are taken from the RC2 
catalog (de Vaucouleurs et al. 1976). Detailed studies show 
that for spiral galaxies the ratio of the optical to radio disk 
scale length is ~ 0.9-1.3 (Hummel 1980; Fitt & Alexander 
1993). For both the normal galaxies (open squares) and the 
starbursts (filled squares), we have computed B m [ n using the 
radio luminosity at 1 .4 GHz. 

The open triangles illustrate the uncertainties in the results 
of Figure [l\at high E g . These points denote the "extreme" 
starbursts of Downes & Solomon (1998) (see their Table 12): 
IC 883 (Arp 193), Mrk 273, and the two individual nuclei of 
Arp 220. For these systems, we use the gas surface densities 
and radial sizes from Downes & Solomon (instead of Kenni- 
cutt 1998), and we combine this data with radio fluxes from 
the 8.44GHz observations of Condon et al. (1991), which 
have sufficient resolution to probe the same spatial scales 
as the observations of Downes & Solomon (in contrast to 
the 1 .4 GHz observations used for the rest of the systems in 
Fig [D- The properties of these starbursts are collected in Ta- 
ble[5] Although the inferred Z? m ; n is somewhat larger at a given 
E g than for the starburst sample in Table [2] (filled squares in 
Fig- ID* the reasonably close correspondence is encouraging. 

In addition to the inferred minimum energy magnetic field, 
Figure ^ a ls° includes several theoretical curves. The dotted 
line shows the equipartition magnetic field (eq. Q) while the 

dashed line shows the scaling B m i n oc E^ 5 derived in 34.11 Fi- 
nally, the solid lines show the scalings B oc E" (with a = 0.9, 
0.8, and 0.7) used in our discussion of the FIR-radio corre- 
lation in 34.31 Although these alternative (non-equipartition) 
scalings are somewhat arbitrary, the particular choice B oc E°- 7 

8 B m in oc ftjj 7 so uncertainties in h ra d introduce only modest uncertainties 
in B mm . 



is physically motivated by setting the magnetic energy den- 
sity equal to the pressure in the ISM produced by star forma- 
tion (P*). Because P+ oc E* (where E* is the star formation 
rate per unit area; see, e.g., Chevalier & Clegg 1985, Thomp- 
son, Quataert, & Murray 2005) and because the Schmidt law 
for star formation is E* oc E„- 4 (Kennicutt 1998), the scaling 
B oc E° 1 follows. 

Our estimates of B m ; n in Figure [2 for normal star-forming 
galaxies are in good agreement with those estimated by other 
authors (e.g., Fitt & Alexander 1993; Niklas 1995; Beck 
2000). When these minimum energy fields can be compared 
with other magnetic field estimates (e.g., from Faraday ro- 
tation) they are generally found to be consistent (e.g., Beck 
2000; Vallee 1995). In addition, the crucial assumption of the 
minimum energy argument, namely that the energy density in 
the magnetic field is roughly equal to the energy density in 
cosmic rays, is confirmed locally in the Galaxy using 7-ray 
observations (Strong et al. 2000). 

In keeping with previous studies, our results in Figure ^ 
show that in normal star-forming galaxies, B m [ n ~ B eq . How- 
ever, the striking result from Figure ^ is that the starburst 
galaxies systematically have B m ; n < B eq . The ratio r/rain = 
(B^ nin /87T)/f is a strongly decreasing function of E^; for a 
typical starburst such as M82, ~ 10~ 2 , while for Arp 220, 
n ■ ~ 10" 4 9 

'/mm lyJ 

If the minimum energy magnetic field strength estimate is 
applicable, the direct implication of Figure ^ is that in star- 
bursts, magnetic fields are dynamically weak compared to 
gravity in the phase of the ISM in which the radio emitting 
electrons reside. Because of the potential importance of this 
conclusion, it is worth critically examining the applicability 
of the minimum energy estimate in the context of luminous 
starbursts (see Beck & Krause 2005 for a more general discus- 
sion). The minimum energy argument rests on the assumption 
that the magnetic energy density is comparable to the cosmic 
ray energy density. Given that the magnetic field and cosmic 
rays could have rather different energy generation and dissi- 
pation mechanisms ("sources and sinks"), the generic appli- 
cability of this assumption is unclear. 

In the context of local spirals, one argument in favor of 
the minimum energy estimate is precisely that it implies a 
total cosmic ray + magnetic energy density in approximate 
equipartition with gravity. If B B m i n , the net cosmic ray + 
magnetic energy density would exceed the gravitational bind- 
ing energy of the gas, and the field and relativistic particles 
would presumably escape from the galactic disk (e.g., Parker 
1965, 1966; Duric 1990). This provides a natural mechanism 
by which the magnetic field strength can adjust to a value of 
~ B m in- However, such an argument does not apply in star- 
bursts given that the minimum energy argument itself implies 

??min<Cl (Fig.Q). 

A more significant worry concerning the applicability of the 
minimum energy argument in starbursts is that strong syn- 

9 Free-free absorption may suppress the observed radio flux density at GHz 
frequencies in dense starbursts (e.g., Condon et al. 1991). However, even in 
our most extreme case, the ULIRG Arp 220, this correction amounts to a 
small increase (~ 10%) in our inferred B m j„. Another uncertainty is the con- 
tribution of free-free emission to the observed radio flux. At GHz frequencies 
the thermal fraction is of order 10% for normal star-forming galaxies (e.g., 
Niklas 1995). For starbursts the thermal fraction is less certain. In any case, 
subtracting the thermal contribution to the radio luminosity would decrease 
the inferred B m j„, strengthening our conclusion that B mm <C B e q in starbursts. 
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chrotron and inverse Compton (IC) losses could lead to rapid 
electron cooling. If the steady state electron energy density 
is significantly lower as a result, then the minimum energy 
magnetic field is an underestimate. 

In the following sections we critically review indepen- 
dent observational and theoretical constraints on the magnetic 
fields in starbursts. We begin by summarizing arguments in 
favor of B ~ B^n and against the hypothesis that magnetic 
fields in starbursts are significantly larger. Several of these ar- 
guments draw heavily on existing results in the literature, but 
are worth reviewing in the present context. Rebuttals to these 
arguments are given in 5@] where we also provide independent 
arguments in favor of B ^> B ra j n . 

3. ARGUMENTS FOR B ~ B min IN STARBURST GALAXIES 

3.1. Zeeman Measurements 

To our knowledge, there are no extragalactic detections 
of Zeeman splitting. In four ULIRGs, however, Killeen et 
al. (1996) derive upper limits of B < 3-5 mG using OH 
masers. These limits are larger than the minimum energy es- 
timate of ~1 mG for ULIRGs, but are a factor of ^5 smaller 
than B eq for Arp 220 in Figure ^ In our own Galaxy OH 
masers trace regions of higher than average surface density 
and magnetic field strength, with B « 2- 10 mG (e.g., Fish et 
al. 2003); it is unlikely that the opposite is true in ULIRGs, in 
which case the Killeen et al. (1996) upper limits argue against 
B ~ B eq . However, a larger sample of Zeeman measurements 
is needed, particularly for galaxies which also have reliable 
gas surface density measurements. 

3.2. Cooling Breaks 

The synchrotron cooling timescale for cosmic ray electrons 
emitting at frequency v is 

^lO^^yr, (4) 

where Z?ioo = B/100fiG. IC losses can be appreciable as well 
and lead to cooling on a timescale 

nc « 5 x 10 5 B\g v-^i U~l_ 9 yr, (5) 

where t/ p h = 10~ 9 £/ p h -9 ergs cm" 3 is the energy density of (pri- 
marily infrared) photons in the starburst. 

The cooling timescale for electrons can be indirectly con- 
strained using the radio spectrum. For continuous injection of 
relativistic electrons, if the cooling time 

r cool = r syn + r IC (6) 

of the electrons emitting at frequency v is shorter than the 
escape time from the galactic disk, the radio spectrum will 
steepen by Aa = 1/2 at higher frequencies, where F v oc v~ a . 
For canonical electron power-law indices of p = 2 expected 
theoretically for strong shocks (e.g., Blandford & Eichler 
1987) and observed in situ in some supernova remnants (e.g., 
Aharonian et al. 2005; Brogan et al. 2005), a = 1/2 in the ab- 
sence of cooling and steepens to a = 1 in the presence of cool- 
ing. The escape time of relativistic particles from the galactic 
disk is uncertain. In local spirals the escape is due to diffusion 
across the magnetic field, while in starbursts it may be much 
more rapid due to advection out of the galaxy with a galactic 
wind. The escape time in the latter case is roughly 

Tesc ~ h/v w - 3 x 10 5 h m Vgo yr, (7) 



where h = 100/zioo pc is the scale height of the galactic disk 
and v w = 500 v'5oo km s" 1 is the speed of the galactic wind (see, 
e.g., Martin 1999 for estimates of the latter). 

The radio spectra of both normal spirals and luminous star- 
bursts are similar with (a) w 0.75 from ss 1 — 10 GHz (e.g., 
Condon et al. 1991; Condon 1992; Niklas et al. 1997). There 
appears to be no strong variation of a with the luminosity 
of the galaxy (though Niklas et al. 1997 note a weak trend 
for starbursts to have somewhat flatter nonthermal spectral in- 
dices). The lack of a systematic variation of the spectral index 
with luminosity, together with the fact that the observed radio 
spectra are significantly flatter than expected for a "cooled" 
electron distribution, suggests that r coo i > r esc , even in lumi- 
nous starbursts. Using equation to estimate T esc , equation 
® then implies field strengths <2x 10~ 4 G, reasonably con- 
sistent with the minimum energy estimates for the starbursts 
in Figure^and inconsistent with B ^> B min . 

3.3. Synchrotron Halos 

The spatial variation of the radio spectrum can provide a 
more stringent constraint on the importance of electron cool- 
ing than the integrated radio spectrum alone. To take a con- 
crete example, observations of M82 by Klein et al. (1988) 
and Seaquist & Odegard (1991) reveal an extended radio halo 
along the minor axis of M82, coincident with the observed 
galactic wind. The radio spectrum at GHz frequencies is 
observed to be roughly constant in the disk of M82, but to 
steepen significantly in the halo a distance h\, lel ± ~ 300 pc 
above the midplane. The usual interpretation of this halo is 
that relativistic electrons generated in the starburst are ad- 
vected out with the galactic wind. The steepening at /ibreak ~ 
300 pc can then be modeled as a consequence of electron 
cooling (e.g., Seaquist & Odegard 1991). The radio emitting 
electrons in the halo must thus have /ztreak ~ t coo [V w . This ap- 
proximate equality can be used to put a limit on the magnetic 
field in the disk of M82. For example, the equipartition field 
of B eq ks 1.6 mG in M82 (see Fig. [3 implies a synchrotron 
cooling time of r syn ~ 10 4 yrs, which, in turn, requires an un- 
physically large v w ~ 30,000 km s" 1 to explain the observed 
scale of the spectral break in M82. Thus a magnetic field 
strength of ~ B eq is strongly disfavored under the interpreta- 
tion that the synchrotron halo in M82 is due to electrons ac- 
celerated in the disk and advected out in a galactic wind. By 
contrast, a field strength of B ~ B m [ n implies v w ~ 500 km s -1 
to explain a spectral break at /ibreak ~ 300 pc — much more 
reasonable in the context of galactic winds. Thus, the spatial 
variation of the radio spectrum in M82 appears to provide in- 
dependent evidence that r coo i > r esc in the galactic disk and 
hence that B ~ B m i n <C B eq . 

In fact, more careful estimates account for both IC cool- 
ing and adiabatic cooling of relativistic electrons as they are 
advected out into the halo. This increases the required wind 
speed by a factor of few to ~ 2000 km s" 1 (e.g., Seaquist & 
Odegard 1991). Although this is substantially smaller than 
the speed of ~ 30, 000 km s" 1 required for B ~ B eq , it is larger 
than the inferred speeds of galactic winds (e.g., Martin 1999) 
and is comparable to the initial velocity of supernova ejecta. 
One cannot, however, rule out a very low density phase of a 
galactic wind with the required velocity ~ 2000 km s -1 . 

We note that NGC 253 (Carilli et al. 1992), NGC 1569 (Is- 
rael & de Bruyn 1988; Lisenfeld et al. 2004), NGC 4631 (Ek- 
ers & Sancisi 1977), NGC 891 (Allen, Sancisi, & Baldwin 
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1978), and NGC 4666 (Dahlem et al. 1997) all exhibit radio 
halos qualitatively similar to those in M82. Unfortunately, 
with the exception of NGC 253, these systems have lower 
gas surface densities than M82 and lie in the part of Figure [2 
where B m ; n ~ B eq . Clearly, a more extensive survey for low 
surface brightness extended radio emission around starbursts 
with Eg ~ 1 g cm" 2 would be very useful in discriminating be- 
tween the minimum energy and equipartition magnetic field 
strengths (see Table|2). 

4. ARGUMENTS FOR B > B min IN STARBURST GALAXIES 

In this section we r eb ut th e arguments for B ~ B m [ n pre- 
sented in Sections f3.2l and l3~3l and provide additional evidence 
that B 3> Bmin- Most of this section centers on showing that 
T C ooi < Tesc is consistent with radio observations of starbursts. 

4. 1 . Understanding B m { n as a Function of Tig 

A magnetic field strength of B ^> B m i n in starbursts would 
imply that the cooling time for relativistic electrons is very 
short compared to the escape time (r coo i <C T esc ). To show 
this, we assume that B = r/ 1//2 Beq, that T esc = 10 6 T eS c.6 yr, and 
we combine equations (0 and (@) to estimate that 

^-O.OIE^V^V^. (8) 

'esc 

For Eg > 0.05 t^!\ rf 1 ^ 1 g cm" 2 , equation (|8} implies that 
Tsyn < Tesc for electrons emitting at GHz frequencies. Includ- 
ing IC cooling decreases this critical surface density by a fac- 
tor of order unity. In addition, for lower surface density galax- 
ies, the electron escape time becomes increasingly determined 
by diffusion across the magnetic field rather than escape in a 
galactic wind. This likely increases r esc substantially and thus 
decreases the critical surface density where r syn ~ r esc . 

The above considerations show that if B ^S> B m \ n then nearly 
all of the starbursts in Figure ^ are in the limit where t coo i < 
T esc . Our observationally inferred results for B m j n (E g ) can 
be understood quantitatively in this limit as follows. In the 
rapidly cooling limit, the electrons radiate all of the energy 
supplied to them by supernova shocks and thus the radio flux 
per unit volume is determined by the supernova rate per unit 
volume: 

vL v /V oc t^/h, (9) 

where E+ is the star formation rate per unit area and V is the 
radio emitting volume. Using the Schmidt Law for star for- 
mation (E* oc Eg 5 ; Kennicutt 1998) and equation Q, we find 

B min oc E 2 / 5 /!" 2 / 7 . (10) 
To estimate the proportionality constant in equation 1 11 PL we 
assume that 10~ 3 £ of each SN's energy is supplied to relativis- 
tic electrons (and is radiated away), in which case equation 
l llOt becomes 

B mm « 7 x l0- 5 6 2 2 /7 e /7 Z 2 g /5 fcloo^GHz G, (ID 
Taking h ~ constant and £ = IP (see 34. 2> . we plot equation 
( II Q in Figure [J (dashed line) and find that it is in reason- 
able agreement with the results for starbursts. This provides 
strong support for the hypothesis that t coo i < r esc in these sys- 
tems. Furthermore, the agreement between the starburst data 
and equation M \\ emphasizes that the inferred minimum en- 
ergy field in starbursts is not a reliable magnetic field mea- 
surement, but rather simply reflects the fact that the radio flux 
is proportional to the star formation rate (eq. |9|). 



4.2. The FIR-Radio Correlation 

Equation (|9} implies a linear correlation between the FIR 
luminosity (Lfir) and the radio luminosity (L ra d) of galaxies. 
Such a correlation is in fact observed over more than four or- 
ders of magnitude in Lfir (van der Kruit 1971, 1973; Helou 
et al. 1985; Condon 1992; Yun, Reddy, & Condon 2PP1). 
Equation is a consequence of the fact that if r coo i < r esc , 
then the radio flux is essentially independent of the magnetic 
field strength and depends only on the rate at which energy 
is supplied to the relativistic electrons. The linearity of the 
FIR-radio correlation then follows naturally from the fact that 
the supernova rate is proportional to the star formation rate 
(e.g., the "calorimeter" theory; Volk 1989; see also Bressan 
et al. 2PP2). Thus as long as r coo i < r esc one can account for 
the observed linear FIR-radio correlation with no fine tuning. 

By contrast, if r coo i > T esc (as argued in 33.21 & 13. 3> then 
significant fine tuning is required to explain the FIR-radio cor- 
relation. Because the synchrotron cooling timescale (eq. |4)) 
varies by a factor of ~ IP 3 for the systems in Figure^(a factor 
of ~ 6P for the starbursts alone), the escape timescale would 
have to vary by exactly the same factor in order to preserve a 
linear FIR-radio correlation. We believe that this level of fine 
tuning argues strongly against t coo i > r esc (and thus indirectly 
against B ~ B min ). 

The magnitude of the radio flux from star-forming galax- 
ies is also fully consistent with the hypothesis that t coo i < 
r esc . To demonstrate this, we note that if each supernova 
supplies an energy of £1P 48 ergs to cosmic ray electrons 
with a spectral index p = 2, and if t coo i <C T esc , then the ra- 
dio luminosity should be related to the FIR luminosity via 
vh v « lp- 5 £L /fl /21n( 7mflx ) « 2.5 x lP" 7 ^, where 
the maximum Lorentz factor of the accelerated electrons. For 
comparison, the 1.4 GHz FIR-radio correlation from Yun et 
al. (2PP1) is i'L l ,«2x 10~ 6 L !R . This implies that £ w 8 
is needed to account for the FIR-radio correlation. That is, 
w P. 8% of each supernova's energy is radiated away as syn- 
chrotron radiation. This is similar to standard estimates for the 
fraction of SN energy supplied to relativistic electrons (see, 
e.g., §2.2.2 of Keshet et al. 2PP3 for a summary), which sup- 
ports the hypothesis that r coo i < r esc . In particular, this argu- 
ment rules out the possibility that t coo i ^> r esc , since in this 
case the amount of supernova energy supplied to cosmic rays 
would have to be prohibitively large (f 3> 1). 

An oft-cited objection to the calorimeter model for the FIR- 
radio correlation is that the nonthermal radio spectra of galax- 
ies should then be significantly steeper tha n are observed be- 
cause of st rong synchrotron cooling (see 33.21 e.g., Condon 
1992). In 34.31 we show that ionization and bremsstrahlung 
losses can systematically flatten the radio spectra of starbursts 
even when r coo i <C T esc , providing much better agreement with 
observations. 

Equation © implies that T syn > r esc for E ? < P.PP5 (P.P5) 
g cm -2 , if r esc f»3 x IP 7 (IP 6 ) yr. The former value of T eS c 
is comparable to that estimated for the Milky Way (Garcia- 
Munoz et al. 1977; Connell 1998). According to the argu- 
ments presented in this section, the calorimeter theory for the 
FIR-radio correlation thus does not apply to low surface den- 
sity normal star-forming galaxies. One might then expect 
a change in the FIR-radio correlation at low E g because a 
smaller fraction of the electron energy is radiated away. Ob- 
servationally, there is only a very mild suppression in the ratio 
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FIG. 2. — The ratio of the minimum energy magnetic energy density to the photon energy density (f/fi min /U p i,) as a function of the gas surface density (Eg). 
The symbols are the same as those in Fig.0 U p i, is calculated using eq. 1121 (see also eq. 1131 ) for all systems except the Galaxy, for which we take F = L/ttR 2 
with L= 1.5 X 10 10 Lq and R = 8.5 kpc (as appropriate for the estimate of B ra j„ at the solar circle). Ub < U p i, is difficult to reconcile with the linearity of the 
FIR-radio correlation and suggests than the minimum energy field strength is an underestimate for luminous starbursts. 



of the radio to FIR luminosities for low-luminosity galaxies 
(Yun et al. 2001). Bell (2003) has shown, however, that this 
is in part because neither the radio nor the FIR emission from 
low luminosity galaxies are good proxies for the total star for- 
mation rate. Because the focus of this paper is on starburst 
galaxies, we defer a more detailed analysis of the connection 
between the FIR-radio correlation in normal spirals and star- 
bursts to a future paper. 

4.2.1. Synchrotron versus Inverse Compton Cooling 

Regardless of the value of T coo i/r esc , the linearity of the 
FIR-radio correlation suggests that the synchrotron cooling 
timescale r syn (eq. |4]) must be shorter than the IC cooling 
timescale t\c (eq. 0). This in turn implies that Ub > Uph- If 
this constraint is not satisfied, any variation in Ub/U p i, would 
imply large changes in the fraction of cosmic ray electron en- 
ergy radiated via synchrotron radiation. A linear FIR-radio 
correlation would then require significant fine tuning. 

In Figure|2]we plot the ratio {Bl^ m /%n)/U p h versus gas sur- 
face density for the sample of galaxies from Figure ^ (Tables 
^|3j. The photon energy density is computed using 

U ph = - = et+c, (12) 
c 



where F is the radiative flux, is the star formation rate 
per unit area, and we take e = 3.8 x 10~ 4 for consistency 
with Kennicutt (1998). Above the dotted line in Figure |3 
UB min /U p h > 1, T syn < tic, and synchrotron dominates IC as 
the primary coolant for the relativistic electrons. Below the 
dotted line r syn > tjc and IC dominates. 

Figure |2 shows that for high surface density starbursts the 
minimum energy argument implies magnetic energy densities 
smaller than the photon energy density. For the most extreme 
case of Arp 220 (and ULIRGs in general), the minimum en- 
ergy estimate implies Us/U p h ~ 0.1. For most starbursts in 
Figure|2l a modest increase of B by a factor of ~ 2 - 3 would 
lead to Ub > U p h- This is plausibly within the uncertainties of 
the minimum energy estimate — and is a rather minor correc- 
tion on the scale of B eq /B m \ n in Figure ^ — but does suggest 
that Bmin is an underestimate. 

Equation H2i may significantly underestimate U p h for 
galaxies with S f > 0.5 g cm" 2 because the ISM is optically 
thick to its own FIR radiation in these systems. In this case, 

U p h = ^— «T /fi eE*e, (13) 
c 

where tjr « KmT, g is the optical depth of the galactic disk to 
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the IR radiation and kjr » 3 — 10 cm 2 g -1 is the temperature- 
dependent dust opacity in the FIR (Thompson, Quataert, & 
Murray 2005). In principle, this factor of tjr could increase 
U p h by a factor of 10- 100 for ULIRGs. If one requires that 
Ub > Up/, in order to satisfy the FIR-radio correlation, the 
necessary field strength would then be 3> B m \ n , and perhaps 
as large as ~ B eq . The applicability of equation (1131 hinges, 
however, on the assumption that the radio emitting electrons 
are co-spatial with optically thick gas at the average density 
of the ISM. If the radio emitting electrons are in a lower-than- 
average density (optically thin) environment, then equation 
O is still valid. 

4.3. Cooling Breaks 

In 33.21 we noted that because starburst radio spectra do 
not show significant evidence for synchrotron cooling, the 
straightforward interpretation is that t coo i > T esc and, hence, 
B ~ B m [ n . One difficulty with this conclusion is provided by 
ULIRGs such as Arp 220, where the IC cooling time alone is 
~ 10 4 yrs. Unless h < 10 pc, the cooling time in ULIRGs 
is shorter than the advection time across the galactic disk, 
for reasonable values of the advection velocity. Given that 
the size of the radio emitting region resolved by the VLA 
in ULIRGs is typically — lOOpc (Condon et al. 1991), one 
might thus expect the nonthermal spectrum to be systemati- 
cally steeper in ULIRGs than in other starbursts. 10 This does 
not however, appear to be the case (Condon et al. 1991; see 
also Fig. 19 of Downes & Solomon 1998 for the case of Arp 
220). The lack of a systematic difference between the radio 
spectral indices of ULIRGs and other systems raises signifi- 
cant questions about the extent to which the radio spectra can 
be easily used to constrain the electron cooling time and thus 
the magnetic field strength. 

We suggest that the radio spectra of starbursts are in fact 
compatible with r coo i <C r esc because of the effect of ioniza- 
tion and bremsstrahlung losses on the low energy part of the 
electron distribution function. The timescale for a relativis- 
tic electron that emits synchrotron radiation at frequency v to 
lose its energy ionizing neutral hydrogen of density n = 10«io 
cm" 3 is 

T lon «10Vo4Hz^OoV- (14) 

Relativistic electrons also lose energy to bremsstrahlung 
emission on ambient neutral nuclei, on a timescale 

T b rem«3xl0 6 «7 ( iyr. (15) 

Setting B = rj^ 2 B eq and n = f(n) = fS g /2hm p we find that 
the ratio of the ionization loss time to the synchrotron cooling 
time is 

— «1.6t/ G Hz/Jioory /2 (16) 

7"syn 

and that the ratio of the bremsstrahlung loss time to the syn- 
chrotron cooling time is 

^^2.4^/z.oor 1 ^ 4 ^ 2 , (17) 

7"syn 

where /1100 = ft/100 pc is the gas scale height. Equations d!6i 
& (I17> show that if the cosmic ray electrons interact with 

10 Although the radio sizes are ~ 100 pc, a small scale height for ULIRGs 
is not ruled out: for a random velocity Sv, the scale height of a galactic disk is 
h rj <5v 2 /27rGE, which implies h 10 pc for the observed values of Sv 100 
km s" 1 and £ « 10 g cm" 2 in Arp 220 (e.g., Scoville et al. 1997). 
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FIG. 3. — Predicted radio spectral index as a function of frequency. Elec- 
trons are injected with ra(-y) ex -y - ?, with p = 2 (solid lines) and p = 2.5 (dashed 
lines), and lose energy through synchrotron radiation, bremsstrahlung, and 
ionization losses. Escape is neglected, which corresponds to the limit 
Tcooi Tesc- The calculations assume T, g = 1 g cm -2 , T] = 1, /jioo = 1. and 
/ = 1,0.5,0.25 (from bottom to top). Ionization losses flatten the radio spec- 
trum at low frequencies. Therefore, even in the limit r coo i <C T eS c, spectral 
indices at GHz frequencies are less than the canonical value of a = 1 (for 
p = 2) appropriate for strong synchrotron cooling. This provides much better 
agreement with the observed radio spectra of starbursts. 



gas at the mean density of the ISM (/ ~ 1) ionization and 
bremsstrahlung losses can significantly modify the distribu- 
tion function for electrons emitting at ~ GHz frequencies. 

Under the action of ionization losses alone, low energy cos- 
mic rays lose energy more rapidly than high energy cosmic 
rays (the bremsstrahlung loss time is independent of electron 
energy). The opposite is true for synchrotron losses. For 
this reason, ionization losses can systematically decrease the 
spectral index from what one would expect for a synchrotron- 
cooled electron distribution (a ~ 1 - 1.25). To demonstrate 
this explicitly, we have calculated the steady state electron 
distribution function for a model in which relativistic elec- 
trons are injected with a spectral index p and lose energy via 
synchrotron radiation, bremsstrahlung losses, and ionization 
losses. Escape is neglected, which corresponds to the limit 
Tcooi "C T esc . Figure |5] shows the resulting synchrotron spec- 
tral slope a as a function of frequency. The solid lines are for 
p = 2 and the dashed lines are for p = 2.5. The calculations 
shown in Figure[3]assume T, g = 1 g cm -2 , 77 = 1, /1100 = 1, and 
/ = 1 , 0.5, and 0.25 (from bottom to top). 

The spectral indices in Figure [3] at GHz frequencies are 
compatible with the observed values even in the presence of 
strong cooling. We find similar results for reasonable varia- 
tions of our model parameters. For example, if p = 2, 77 = 0.1, 
and / = 0.25, we find a w 0.6 at 1 GHz from E g w 0.1 - 10 
g cm -2 and a « 0.7-0.8 at 5 GHz over the same range of 
T, g , Our calculations predict that the non-thermal spectral 
index should steepen to a ~ 1 at sufficiently high frequen- 
cies, but this may be difficult to detect because of free-free 
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emission. A second prediction of this model is that the ra- 
dio spectral index should be larger at high frequencies than it 
is at low frequencies, with a transition at ~ GHz frequencies 
where the ionization, bremsstrahlung, and synchrotron loss 
times are comparable (see eqns. 1161 & 1171 ). For example, 
in the model from Figure [5] with p = 2 and / = 0.5 (the mid- 
dle solid line), the effective spectral index from 1-10 GHz 
is a = 0.8, whereas from 0.1-1 GHz it is a = 0.55. Using 
the radio data in Table|2 as well as 57.5, 151, 408, 750, and 
4850 MHz data available in the NASA Extragalactic Database 
(from Hales et al. 1993; Ficarra et al. 1985; Israel & Ma- 
honey 1990; Becker et al. 1991; McCutcheon 1973; Waldram 
et al. 1996; Heeschen & Wade 1964; Large et al. 1981; Grif- 
fith et al. 1995; Kuehr et al. 1981), and additional data from 
Oly & Israel (1993) (for NGC 3504), Irwin & Saikia (2003) 
(for NGC 3079), and Hummel et al. (1991) (for NGC 891), 
we find that for the sample of starburst galaxies considered in 
this paper the spectral index below 1 GHz is indeed typically 
flatter than that above 1 GHz by Aa w 0.3 (although there are 
some exceptions and scatter). 

Although encouraging, a larger and more homogeneous 
sample of data is clearly needed to test our predictions of 
spectral index variations as a function of frequency. In addi- 
tion, low frequency radio emission may be strongly modified 
by free-free absorption in many starbursts. Moreover, elec- 
trons emitting at ~ 100 MHz in a ~ mG field have energies of 
^100 MeV. At these energies secondary electrons produced 
by charged pion decay are important and may significantly 
modify the electron spectrum as well. 

Because the synchrotron, ionization, and bremsstrahlung 
loss rates all scale differently with the gas density and mag- 
netic field strength, there is no guarantee that the linearity of 
the FIR-radio correlation is preserved under the influence of 
ionization and bremsstrahlung losses. To assess this, Figure 
|4]shows model calculations of LpiR/Zrad as a function of E g . 
The synchrotron flux at 1 GHz is calculated assuming that 
electrons are injected with a p = 2 spectrum, and lose energy 
via synchrotron, ionization, and bremmstrahlung losses (as 
before, escape is neglected). The FIR luminosity is assumed 
to be proportional to the rate at which energy is supplied to 
relativistic electrons. The models shown in Figure |4] assume 
that / = 0.5 (n = 0.5 (n)), h m = 1, and that BaE" with a = 
1,0.9,0.8, and 0.7 (bottom to top; see Fig. [Hforthe B(T, g ) 
curves). All of the assumed models for B(E ? ) have B ^> B m \ n 
for luminous starbursts. 

Figure|4]demonstrates that the predicted FIR-radio correla- 
tion is linear to better than a factor of 2 over 3 orders of mag- 
nitude in Eg. The models with a = 0.7, 0.8, and 0.9 provide 
particularly good agreement with the observed linearity of the 
FIR-radio correlation. Because galaxies of a given E g can 
have a range of FIR luminosities (depending on galaxy size), 
the systematic trends in Figure[4]may show up in part as scat- 
ter in the canonical plots of L^i^/L^d vs. Lfir; the < factor of 
2 variation in Figure |4]is consistent with the observed scatter 
in the FIR-radio correlation. The approximate constancy of 
LpiR/Lrad in Figure 0] arises because if / ~ 1 and B » B m in, 
then ^nAik^ + ibrem) -1 is relatively constant over most of the 
parameter space relevant to observed starbursts. 

In ^4.2l we estimated that » 0.8% of the energy of each su- 
pernova must be supplied to relativistic electrons to account 
for the magnitude of the w 1 .4 GHz radio flux from starburst 
galaxies, taking into account synchrotron radiation alone and 
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FIG. 4. — Predicted ratio of the FIR to 1 GHz synchrotron luminosity as 
a function of gas surface density S g . Electrons are injected with a p = 2 
spectrum and lose energy through synchrotron radiation and ionization and 
bremsstrahlung losses; escape is neglected. The FIR luminosity is assumed 
to be proportional to the total rate at which energy is supplied to relativistic 
electrons. The calculations assume that / = 0.5 (n = 0.5 (n)) and that B oc S° 
with a = 1, 0.9, 0.8, and 0.7 (bottom to top) and B = B cq at S g = 0.01 g 
cm -2 (see Fig.^for these scalings). These results demonstrate that although 
ionization and bremsstrahlung losses modify the radio spectra of starbursts 
(Fig. |3}, they nonetheless maintain a nearly linear FIR-radio correlation, in 
agreement with observations. 



assuming electrons are injected with a p = 2 spectrum. If ion- 
ization and bremsstrahlung losses are important in starbursts, 
as argued in this section, then the required energy injection 
rate increases to w 1.6% of the supernova energy. In addition, 
if the electron injection spectrum is p = 2.2, then the required 
energy injection increases by an additional factor ps 2. 

In contrast to the results of Figure |4] in a model in which 
/ = 0.5 and B oc E^ , which corresponds to the observed scal- 
ing of B m j n (E f ) in Figure n the ratio of the FIR to radio lu- 
minosities increases by a factor of w 10 from E g = 0.01 g 
cm" 2 to E ? = 10 g cm" 2 because of the increasing suppres- 
sion of the radio flux by bremsstrahlung losses. Thus field 
strengths w Z? m ; n are only consistent with the observed radio 
emission from starbursts if cosmic-ray electrons interact with 
gas at much less than the mean density of the ISM (/ <C 1). 

To summarize this subsection, if B 3> B m ; n and if cosmic- 
ray electrons interact with gas at approximately the mean den- 
sity of the ISM, ionization and bremsstrahlung losses will sys- 
tematically flatten the nonthermal spectra of starburst galaxies 
(Fig. |3j andyet maintain the linearity of the FIR-radio corre- 
lation (Fig.|4). This eliminates one of the strongest arguments 
against t coo i <C r esc and argues for B » B m \ n in starbursts. 

4.4. Pion Production 

An independent test of the density of the ISM in which 
the cosmic rays reside can be provided by gamma-ray emis- 
sion from neutral pion decay. The timescale for cosmic ray 
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protons to lose their energy to pion production is f p i on ~ 
10 8 («/lcm" 3 ) -1 yr. If n ~ (n), this is significantly shorter 
than r esc in starbursts (see eq. Q) an d so pi° n production will 
lead to a substantial 7-ray luminosity. Assuming that 10% of 
the energy of each supernova is supplied to relativistic pro- 
tons with a spectral index p = 2, and that w 1 /3 of this en- 
ergy goes into 7-rays rather than neutrinos or pairs, we esti- 
mate a 7-ray luminosity from neutral pion decay in starbursts 
of vL v k 3 x 10" 4 L/R/ln(7„, flA ) w 2 x 10~ 5 L /R , where 7„, flA is 
the maximum Lorentz factor of the accelerated protons and 
the estimate is only applicable above w 100 MeV (see Tor- 
res 2004 for more detailed calculations). This predicted flux 
is just below the EGRET upper limits for systems such as 
M82 and Arp 220 (Cillis et al. 2005), but should be detectable 
with GLAST. Because the timescales for electron ionization 
and bremsstrahlung losses and proton pion losses are quite 
similar, detection of gamma-ray emission at this level would 
also confirm the importance of ionization and bremsstrahlung 
losses for the observed radio spectra of starbursts. 

Relativistic protons lose « 2/3 of their energy to the pro- 
duction of charged pions, which subsequently decay into high 
energy neutrinos and electron/positron pairs; the pairs receive 
« 1/4 of the pion energy. Thus if ~ 10% of the energy of each 
supernova is supplied to relativistic protons, and if f p ; on -C t esc , 
then ~ 1.5% of the energy of each supernova goes into sec- 
ondary electrons and positrons. This is in good agr eeme nt 
with the energy injection rate estimated in q4.2l and 34.31 to 
account for the radio flux from starbursts. This suggests that 
secondary electrons may dominate the radio emission in star- 
bursts (see also Rengarajan 2005). If correct, it remains to 
be understood why normal spirals and starbursts lie on ap- 
proximately the same FIR-radio correlation even though sec- 
ondaries do not appear to be important for electrons emitting 
at ~ GHz freque ncie s in normal spirals (Strong et al. 2004). 
As discussed in 34.2l variations in T syn /r esc might also be ex- 
pected to modify the FIR-radio correlation for normal spirals 
relative to starbursts. 

4.5. Inverse Compton Emission 

IC upscattering of infrared photons by cosmic-ray elec- 
trons can provide an independent constraint on Ub/U p i, be- 
cause the ratio of the IC power to the radio power is given by 
Ljc/L m d ~ Uph/Ug- If we use equation dl2> to estimate U p h, 
then B ~ B eq in starbursts would lead to negligible IC power 
because it would imply Ub /U p i, 3> 1 ■ However, if the infrared 
optical depth is larger than unity, equation Jl 31 may instead 
be applicable, in which case B ~ Z? eq implies Ub/U p i, ~ 1. 

IC emission has long been argued to contribute significantly 
to the X-ray emission in a variety of systems, including M82 
(e.g., Hargrave 1974; Moran & Lehnert 1997) and NGC 3256 
(e.g., Moran, Lehnert, & Helfand 1999). In both of these ex- 
amples, however, much of the hard X-ray emission has been 
resolved into point sources with Chandra, although IC emis- 
sion may still contribute to the hard X-rays from the very cen- 
tral nucleus of the starburst (e.g., Lira et al. 2002; Strickland 
et al. 2002; Strickland & Heckman, private communication). 
More detailed tests of this hypothesis would be very worth- 
while given the possibility of providing an independent con- 
straint on the magnetic field strength in starbursts. 

4.6. Synchrotron Halos 



As an alternative to the standard interpretation (" 33. 3> that 
extended synchrotron halos are due to relativistic electrons 
generated in the disk and advected into the halo (before cool- 
ing), we suggest that the halos are instead due to electrons ac- 
celerated in situ in the galactic wind, by either direct electron 
acceleration in shocks or via charged pion production. We 
note that in several simulations of galactic winds (e.g., Strick- 
land & Stevens 2000), most of the X-ray emission is due to in 
situ shocks, which are also likely to be sites of particle accel- 
eration. Since the total kinetic energy flux in a galactic wind 
is comparable to the energy produced by supernovae in the 
starburst, only a small fraction of the wind energy must be 
converted into relativistic particles in situ to account for the 
extended synchrotron emission observed. In this interpreta- 
tion, which clearly requires additional work, B 3> B m j n would 
be fully consistent with extended synchrotron emission. Ob- 
served spectral differences between the galactic disk and halo 
in M82 and other systems may arise from spatial variations in 

7"syn/7esc? tlC/f&SCt 7brem/7syn» and T[ on fT S y n . 

5. DISCUSSION & CONCLUSIONS 

Based on the analysis of the preceding sections, it is worth- 
while contrasting the following two scenarios for the radio 
emission from starbursts: 

B ~ fi m i n : For the minimum energy magnetic field estimate to 
be applicable, the electrons must escape rapidly in a galactic 
wind with t coo i > r esc . The rapid escape implies that there is 
no cooling break at GHz frequencies, which is consistent with 
the observed nonthermal spectra of a w 0.75. Rapid escape 
is plausible given that hot outflowing galactic winds are ob- 
served from starbursts and that the cosmic rays are likely gen- 
erated in the same supernova shocks that generate the wind. 
This interpretation of magnetic fields in starbursts — the stan- 
dard one — has two major drawbacks: (1) it is difficult to 
see how r coo i > r esc is applicable in ULIRGs where cooling 
times from IC alone are ~ 10 4 yr and (2) it is very difficult 
to account for the FIR-radio correlation with r coo i > r esc be- 
cause variations in r coo i or r esc should lead to variations in the 
fraction of the electron energy radiated away via synchrotron 
radiation. Either several coincidences or a complex feedback 
loop are then required to account for the linearity of the FIR- 
radio correlation. The latter would be somewhat surprising 
given that B m ; n itself implies that magnetic fields and cosmic 
rays are dynamically weak compared to gravity (Fig.^. 

B 3> B m i n : An alternative possibility is that magnetic fields 
in starbursts are much stronger than is suggested by the min- 
imum energy estimate. Rapid cooling of relativistic electrons 
in starbursts will invalidate the minimum energy estimate if 
Tcooi <C T esc . Magnetic fields in starbursts could then in prin- 
ciple be as large as ~ B eq , which is ~ 10 times larger than 
B m in for typical starbursts (Fig.Q. This model naturally ac- 
counts for the linearity of the FIR-radio correlation because 
in the limit r coo i < r esc the radio flux is nearly independent 
of the magnetic field strength (e.g., Volk 1989). In addition, 
the observed radio luminosities from star-forming galaxies are 
comparable to the estimated rate at which en ergy is supplied 
to relativistic electrons in supernova shocks ( 34.2> . This sup- 
ports the hypothesis that r coo i < T esc . Independent support for 
this hypothesis is provided by the fact that the observed scal- 
2/5 

ing Z? m i„ oc Eg for the minimum energy field in starbursts 
(Fig. ^ follows directly from the Schmidt Law for star for- 
mation when r coo i < r esc (eq. 1101 ). 
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A standard objection to r coo i < r esc is that the typical non- 
thermal spectral indices of starbursts {a « 0.75) do not show 
the expected steepening due to strong cooling (a sw 1). We 
have argued, however, that if B ^> fi m j n and if cosmic ray elec- 
trons interact with gas at approximately the mean density of 
the ISM, then ionization and bremsstrahlung losses flatten the 
radio spectra of starbursts at ~ GHz frequencies and reconcile 
Tcooi S T esc with the observed spectra (Fig. [3}. An important 
part of this argument is the realization that the ionization and 
bremsstrahlung loss times are similar to the synchrotron cool- 
ing time in all starburst galaxies for cosmic ray electrons emit- 
ting at GHz frequencies. Thus, ionization and bremsstrahlung 
losses can modify the nonthermal spectra of starburst galax- 
ies and yet maintain the linearity of the FIR-radio correlation 

(Eg. a. 

This interpretation of magnetic fields in starbursts — that 
B ^> B mm — has two potential drawbacks: (1) Given that 
galactic winds efficiently remove mass and metals from galax- 
ies, it is unclear whether the cosmic rays actually interact 
with the bulk of the ISM in starbursts, which is required for 
ionization losses to be significant. (2) The nonthermal syn- 
chrotron halos observed in several systems are typically in- 
terpreted as cosmic-ray electrons advected out with a galactic 
wind, with synchrotron cooling only becoming important in 
the halo (e.g., Seaquist & Odegard 199 1). If this interpreta- 
tion is correct, it requires B ~ B m ; n ( ^3.3> . We have suggested, 
however, that extended synchrotron emission may be due to 
particles accelerated in situ in galactic winds ( 34.6> . in which 
case the observed halos cannot be readily used to constrain 
the magnetic field strength in galactic disks. 

We believe that the arguments presented in this paper 
strongly favor r coo i < r esc for the relativistic electrons in star- 
bursts. In this limit, the minimum energy argument under- 
estimates the true magnetic field strength. The observed ra- 
dio flux is, however, nearly independent of the magnetic field 
and depends primarily on the rate at which energy is supplied 
to relativistic electrons. Thus direct constraints on the true 
field strength are difficult to come by. By themselves, these 
arguments do not imply that B ~ B eq , only that r coo i < r esc 
and B > B m [ n . Our strongest argument that magnetic fields 
are likely to be ^> fi m j n derives from the ratio of the ioniza- 
tion and bremsstrahlung loss times to the synchrotron cooling 
time: if cosmic ray electrons interact with gas having roughly 
the mean density of the ISM, magnetic fields must be ^> B m \ n 
in starbursts in order for ionization and bremsstrahlung losses 
to not completely dominate over synchrotron losses (which 
would be inconsistent with the FIR-radio correlation and 
would require that an unreasonably large fraction of the su- 
pernova energy be supplied to cosmic ray electrons in orde r 
to account for the observed radio fluxes from starbursts; 34. 3> . 

Our conjecture is thus that the magnetic field is signifi- 
cantly larger than ~ B m i n (the canonical estimate) in many 
starbursts. Confirming this prediction would have significant 
implications for understanding the physics of star formation 
in starbursts, the effects of magnetic fields on the dynamics 
of galactic winds, and the role of magnetic stresses in driv- 



ing gas to smaller radii to fuel a central active galactic nu- 
cleus. Towards this end, we briefly summarize several ob- 
servations that should help test the predictions of our model: 
(1) If ionization and bremsstrahlung losses are important, the 
radio spectra of starbursts shoul d be flatter below ~ 1 GHz 
than they are above ~ 1 GHz ($0]and Fig. 0. (2) If cos- 
mic rays interact with gas at about the mean density of the 
ISM, as is required for ionization and bremsstrahlung losses 
to be important, there should be an appreciable ga mm a-ray 
flux from starbursts due to neutral pion production ( 84.41 . We 
estimate vL v ks 2 x 10~ 5 L/r above ks 100 MeV. This predic- 
tion is testable with GLAST. (3) Zeeman measurements of 
starbursts can directly probe the magnetic field strength in the 
dense phase of the ISM, and distinguish between B ~ B m in and 
B ^> B m i n . Existing upper limits in 4 ULIRGs from Killeen 
al . (19 96) are above B m i n , but below the equipartition field 
( 33. 1> . (4) If B 3> Bmin, IC emission is unlikely t o co ntribute 
significantly to the X-ray emission in starbursts ( 34.5> . 

It is interesting to compare our results on magnetic fields in 
starbursts with the inferred correlation between magnetic field 
strength, column density, and mass density in Galactic molec- 
ular clouds (e.g., Troland & Heiles 1986; Crutcher 1999). 
Using Zeeman splitting to determine the magnetic field 
strength directly, Crutcher (1999) finds that over the range 
0.01 gem" 2 < S H2 < 2gcm" 2 < and 2.7 < log 10 [« H2 (cm" 3 )] < 

1 12 

6.8, that B oc n H and B oc Eh,. These findings imply both 

that that the Alfven speed Va ~ constant ~ 10 km s" 1 and that 
equipartition obtains. Observations of polarization in H2O 
masers in dense star-forming regions provide evidence that 
these trends extend to yet higher densities (10 s — 10 10 cm -3 ; 
Sarma et al. 2002; Vlemmings et al. 2005). These results sug- 
gest that the magnetic field is always comparable to the total 
pressure in Galactic molecular clouds, with B w B eq . 

The sample of Galactic molecular clouds reviewed by 
Crutcher (1999) covers precisely the same surface densities 
represented by starbursts in our Figure^ It is thus particularly 
striking that the minimum energy estimate implies B <C B sq in 
starbursts, while direct Zeeman measurements of local molec- 
ular clouds with the same S s imply B ~ B eq . Our suggestion 
is that this difference arises in part from the inapplicability of 
the minimum energy estimate in starbursts. 
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TABLE 1 

Properties of Normal Star-Forming Galaxies 



Object Name 



D b Iog 10 [E- 



log 10 [S, 



log^l/^dp Angular Diameter Radius fi m ; n B^ 1 Radio Data 



(Mpc) (gcnT 1 ) (Mq yr~ j kpc ~) (mJy) (WHz H ) 



(") 



(kpc) (fiG) (pG) References 



NGC224 ■■■ 0.9 -3.00 

NGC598 •■■ 0.9 -2.65 

NGC628 0.0022 9.3 -2.75 

NGC772 0.0082 34.6 -2.74 

NGC925 •■■ 9.3 -2.77 

NGC 1569 • ■ ■ 2.2 -2.35 

NGC2336 0.0074 31.1 -2.77 

NGC 2403 ••■ 3.3 -2.80 

NGC 2841 0.0021 9.0 -2.71 

NGC 2903 ••■ 6.1 -2.82 

NGC 2976 ••■ 3.4 -2.70 

NGC 3031 ••■ 3.3 -2.83 

NGC 3310 0.0033 13.9 -2.54 

NGC 3338 0.0043 18.3 -2.87 

NGC 3368 0.0030 12.6 -2.75 

NGC 3486 0.0023 9.6 -2.80 

NGC 3521 0.0027 11.3 -2.46 

NGC 3631 0.0039 16.3 -2.52 

NGC 3675 0.0026 10.8 -2.69 

NGC 3726 0.0028 12.0 -2.62 

NGC 3893 0.0032 13.7 -2.62 

NGC 3938 0.0027 11.4 -2.48 

NGC 4178 0.0013 5.3 -2.55 

NGC 4254 0.0080 33.9 -2.29 

NGC 4258 0.0015 6.3 -3.09 

NGC 4294 0.0012 5.0 -2.66 

NGC 4303 0.0052 22.1 -2.47 

NGC 4321 0.0053 22.4 -2.54 

NGC 4394 0.0031 13.0 -3.04 

NGC 4402 0.0008 3.3 -2.60 

NGC 4501 0.0076 32.1 -2.59 

NGC 4519 0.0041 17.2 -2.69 

NGC 4535 0.0065 27.6 -2.67 

NGC 4548 0.0016 6.8 -2.99 

NGC 4569 ••• 18.0 -3.07 

NGC 4571 ••• 15.0 -2.85 

NGC 4579 0.0051 21.4 -2.87 

NGC 4647 0.0047 19.9 -2.64 

NGC 4651 0.0028 12.0 -2.62 

NGC 4654 0.0035 14.6 -2.58 

NGC 4689 0.0054 22.8 -2.74 

NGC 4698 0.0033 14.1 -3.43 

NGC 4713 0.0022 9.2 -2.64 

NGC 4736 0.0010 4.4 -3.03 

NGC 4826 0.0014 5.8 -3.01 

NGC 5033 0.0029 12.4 -2.75 

NGC 5055 0.0017 7.1 -2.51 

NGC 5194 0.0015 6.5 -2.21 

NGC 5236 0.0017 7.3 -1.98 

NGC 5457 ••• 8.0 -2.59 

NGC 6217 0.0045 19.2 -2.39 

NGC 6503 ••• 5.2 -2.79 

NGC 6643 0.0049 20.9 -2.57 

NGC 6946 ••• 3.0 -2.38 

NGC 7331 0.0027 11.6 -2.60 



-3.13 
-2.47 
-2.18 
-2.84 
-2.44 
-0.80 
-1.92 
-2.15 
-2.99 
-2.31 
-1.66 
-2.50 
-1.14 
-2.56 
-2.55 
-2.46 
-1.91 
-1.73 
-2.01 
-2.28 
-1.96 
-2.11 
-2.27 
-1.70 
-2.36 
-1.87 
-1.74 
-2.07 
-2.88 
-2.80 
-2.21 
-1.98 
-2.38 
-2.52 
-2.78 
-2.56 
-2.32 
-2.22 
-1.98 
-2.06 
-2.38 
-3.55 
-1.53 
-2.22 
-2.47 
-2.64 
-2.32 
-1.78 
-1.41 
-2.46 
-1.91 
-2.08 
-1.81 
-1.88 
-2.33 



8400 
3300 
180 
71.4 
46.0 
411 
17.7 
330 
83.8 
407 
50.8 
380 
383 
28.1 
30.0 
59.4 
357 
80.8 
43.7 
27.4 
134 
61.7 
23.3 
422 
790 
24.9 
416 
340 
0.7 
52.2 
278 
9.1 
64.5 
4.6 
83.4 
4.2 
103 
35.0 
24.0 
117 
9.6 
0.6 
37.8 
254 
103 
178 
390 
1490 
2445 
750 
78.4 
36.5 
92.2 
1395 
373 



20.89 
20.51 
21.27 
22.01 
20.68 
20.40 
21.31 
20.64 
20.91 
21.26 
19.85 
20.70 
21.95 
21.06 
20.76 
20.82 
21.74 
21.41 
20.79 
20.68 
21.48 
20.99 
19.90 
22.77 
21.58 
19.88 
22.39 
22.31 
19.15 
19.85 
22.54 
20.51 
21.77 
19.42 
21.51 
20.06 
21.76 
21.23 
20.62 
21.48 
20.78 
19.16 
20.59 
20.78 
20.61 
21.52 
21.38 
21.88 
22.19 
21.77 
21.54 
20.08 
21.69 
21.17 
21.78 



9912 
3354 
612 
432 
534 
150 
414 
948 
408 
642 
294 
1332 
210 
354 
390 
426 
486 
276 
354 
324 
234 
318 
240 
312 
906 
150 
354 
408 
234 
186 
360 
186 
378 
306 
474 
222 
306 
204 
210 
258 
234 
222 
156 
630 
480 
546 
660 
600 
660 
1614 
180 
294 
204 
642 
510 



20.9 
7.3 
13.8 
36.3 
12.0 
0.82 
31.2 
7.6 
8.9 
9.5 
2.4 
10.6 
7.1 
15.7 
12.0 
9.9 
13.4 
10.9 
9.3 
9.4 
7.8 
8.8 
3.1 
25.7 
13.9 
1.8 
19.0 
22.1 
7.4 

I. 5 

28.0 
7.8 

25.3 
5.1 

20.7 
8.1 

15.9 
9.9 
6.1 
9.1 

12.9 
7.6 
3.5 
6.8 
6.7 

16.4 

II. 4 
9.4 

I 1.6 
31.4 
8.4 
3.7 
10.3 
4.6 
14.3 



3.3 
4.6 
5.3 
5.0 
3.9 
15.1 
3.4 
4.9 
5.4 
6.5 
5.6 
4.2 
12.2 
4.3 
4.1 
4.8 
7.4 
6.7 
4.9 
4.5 
8.5 
5.7 
5.1 
10.0 
6.5 
6.8 
9.2 
8.1 
1.9 
7.4 
8.2 
4.5 
5.2 
2.8 
4.9 
3.2 
6.7 
6.2 
5.5 
7.7 
4.0 
1.9 
7.4 
5.8 
5.2 
5.7 
6.4 
9.8 
10.7 
4.6 
8.4 
5.1 
8.2 
9.3 
7.3 



2.3 
5.1 
4.1 
4.2 
3.9 
10.2 
3.9 
3.6 
4.5 
3.5 
4.6 
3.4 
6.6 
3.1 
4.1 
3.6 
7.9 
6.9 
4.7 
5.5 
5.5 
7.6 
6.4 
1 1.7 
1.9 
5.0 
7.7 
6.6 
2.1 
5.7 
5.9 
4.7 
4.9 
2.3 
1.9 
3.2 
3.1 
5.2 
5.5 
6.0 
4.2 
0.8 
5.2 
2.1 
2.2 
4.1 
7.1 
14.1 
23.9 
5.9 
9.3 
3.7 
6.2 
9.5 
5.7 



REFERENCES. — (1) Condon (1987); (2) Condon, Yin, & Burstein (1987); (3) Condon et al. (1990). 

"Redshifts from NED. b Distance. In most cases, distances are computed from z with H = 71 km s _I Mpc -1 . In cases where redshifts are not a reliable distance 
indicator, D is taken from the literature. c Gas surface density. d Data from Kennicutt (1998). c Star formation rate per unit area. f Radio flux density listed at 
1.4 and 1.49 GHz. E Radio luminosity computed using L la( j = 4nD 2 Si ,4ghz- h B m i n calculated using equation Q, with fQHz = 1.49, V = 27rK 2 /2 la( j, and assuming 



Kad = 500 pc. 'B e q calculated using equation |3j. 
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TABLE 2 

Properties of Starburst Galaxies 



Object Name 




£>» 


iog 10 [s s r- d 


l°g 10 [5X] c - d 


Sl.4GHz f 


logiot^ad] 8 


Angular Diameter d 


Radius 


d h 
"min 




Radio Data 






(Mpc) 


(gem -2 ) 


(M yt' 1 kpc- 2 ) 


(mJy) 


(WHz 4 ) 


(") 


(kpc) 


(mG) 


(mG) 


References 


NGC 253 




3.5 


-0.33 


1.24 


5594 


21.93 


24 


0.21 


0.14 


1.1 


1 


NGC 520 


0.0074 


31.2 


0.13 


1.32 


167 


22.29 


5 


0.38 


0.13 


3.1 


2 


NGC 660 


0.0028 


12.0 


-1.08 


0.06 


387 


21.83 


31 


0.90 


0.06 


0.19 


2 


NGC 828 


0.0179 


75.7 


-0.02 


1.10 


108 


22.87 


5 


0.92 


0.11 


2.2 


3 


NGC 891 


0.0018 


7.4 


-1.07 


-0.58 


701 


21.67 


35 


0.63 


0.06 


0.19 


3 


NGC 1097 


0.0043 


18.0 


-1.01 


-0.20 


415 


22.21 


35 


1.52 


0.05 


0.22 


3 


NGC 1614 


0.0159 


67.3 


0.04 


1.79 


123 


22.83 


4 


0.65 


0.13 


2.5 


2 


NGC 1808 


0.0034 


14.2 


-1.03 


0.08 


519 


22.10 


30 


1 .03 


0.06 


0.21 


1 


NGC 2146 


0.0030 


12.6 


-0.85 


0.84 


1087 


22.32 


17 


0.52 


0.11 


0.32 


1 


NGC 2623 


0.0185 


78.0 


-0.81 


1.00 


98.5 


22.86 


8 


1.51 


0.08 


0.35 


2 


NGC 3034 (M82) 




3.3 


-0.16 


1.48 


7657 


22.00 


29 


0.23 


0.14 


1.6 


2 


NGC 3079 


0.0038 


15.9 


0.57 


1.63 


849 


22.41 


5 


0.19 


0.20 


8.5 


2 


NGC 3351 


0.0026 


11.0 


-0.85 


0.24 


47.8 


20.84 


14 


0.37 


0.05 


0.32 


2 


NGC 3504 


0.0051 


21.7 


-0.78 


0.11 


265 


22.18 


16 


0.84 


0.08 


0.38 


2 


NGC 3627 


0.0024 


10.2 


-1.40 


-0.77 


458 


21.76 


39 


0.97 


0.05 


0.09 


2 


NGC 3690 


0.0100 


42.2 


-1.40 


-0.10 


658 


23.15 


24 


2.46 


0.08 


0.09 


2 


NGC 6240 


0.0243 


102.7 


0.43 


1.87 


392 


23.70 


3 


0.75 


0.22 


6.2 


3 


NGC 7552 


0.0053 


22.3 


-1.30 


0.16 


276 


22.22 


22 


1.19 


0.06 


0.1 1 


1 


IC 342 


0.0010 


4.4 


-1.62 


-0.41 


2250 


21.71 


67 


0.71 


0.06 


0.05 


1 


IC 694 


0.0104 


43.8 


0.42 


2.40 


153 


22.55 


2.6 


0.28 


0.18 


6.0 


4 


IC 883 (Arp 193) 


0.0233 


98.5 


0.27 


1.54 


101 


23.07 


3.4 


0.81 


0.14 


4.26 


2 


IC 1623 


0.0204 


86.1 


0.13 


1.67 


221 


23.30 


3.6 


0.75 


0.17 


3.1 


2 


Maffei 2 




3.0 


-1.22 


-0.27 


286 


20.48 


40 


0.29 


0.05 


0.14 


5 


Arp 55 


0.0393 


165.9 


-0.95 


0.32 


33 


23.04 


8 


3.22 


0.06 


0.26 


6 


Arp 220 


0.0181 


76.6 


1.08 


2.98 


324 


23.36 


2 


0.37 


0.26 


27.5 


2 


IRAS 10173+0828 


0.0480 


202.8 


-1.27 


0.48 


9.9 


22.69 


7 


3.44 


0.05 


0.12 


6 


IRAS 17208-0014 


0.0425 


179.4 


0.41 


2.01 


73 


23.45 


3 


1 .30 


0.14 


5.9 


6 


VII Zw 31 


0.0537 


226.8 


-0.57 


0.82 


41.6 


23.41 


5 


2.75 


0.09 


0.62 


3 


Zw 049.057 


. 0.0332 


140.1 


0.22 


1.77 


46.4 


23.04 


3 


1.02 


0.12 


3.8 


2 



REFERENCES. — (1) Condon (1987); (2) Condon et al. (1990); (3) Condon et al. (1996); (4) Condon & Broderick (1988); (5) Condon et al. (1998); (6) Yun, 
Reddy, & Condon (2001). 

"Redshifts from NED. ^Distance. In most cases, distances computed with H = 71 km s Mpc~' . In some cases, distances are taken from the literature instead 
of using the Hubble Law. c Gas surface density. d Data taken from Kennicutt (1998). e Star formation rate per unit area. f Radio flux density listed at 1.40, 1.425, 
or 1.49 GHz. g Radio luminosity computed using L raa = 4nD 2 S\AGHz- h Bmin calculated using equation Q, with vghz = 1.49, V = 2nR 2 h m ^, and assuming 
''rad = 100 pc. 'B e q calculated using equation |3j. 



TABLE 3 

Properties of "Extreme" Starburst Galaxies 



Object Name 




D b 




log l0 [£„] d 


Sg.44GHz 


log 10 [£md]' 


Angular Diameter L 


Radius 


Bmin g 




Radio Data 






(Mpc) 


(gem -2 ) 


(M yr -1 kpc -2 ) 


(mJy) 


(WHz 4 ) 


(") 


(kpc) 


(mG) 


(mG) 


References 


IC 883 (Arp 193) 


0.0233 


98.5 


0.25 


2.68 


34.9 


22.61 


0.69 


0.16 


0.33 


4.06 


1 


Mrk 273 


0.0378 


159.6 


0.66 


3.36 


43.5 


23.13 


0.35 


0.14 


0.51 


10.45 


1 


Arp 220 (west) 


0.0181 


76.6 


0.94 


3.55 


72.2 


22.71 


0.40 


0.07 


0.55 


19.91 


1 


Arp 220 (east) 


0.0181 


76.6 


0.78 


2.96 


60.7 


22.63 


0.66 


0.12 


0.39 


13.77 


1 



References. — (1) Condon et al. (1991). 

"Redshifts from NED. b Distances computed with H = 71 km s _1 Mpc -1 . c Gas surface density from Downes & Solomon (1998), Tables 3, 4, and 12. Values 
adjusted for assumed Hq . d Star formation rate per unit area computed using E* =L m (< R)/(ttR 2 <ec 2 ) with e = 3.8 X 10" 4 (see Downes & Solomon 1998 Table 12). 
c Radio flux density at 8.44 GHz. f Radio luminosity computed using L rat | = 4nD 2 Sg ,44GHz- E #min calculated using equation Q, with ^GHz = 8.44, V = 2nR 2 h iai j, 
and assuming /i la( j = 100 pc. h B C q calculated using equation J3}. 



